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A strong suppression of non-photonic electrons in Au+Au collisions is observed at Rmc. 
This is in disagreement with the expected dominance of the energy loss via gluon radi- 
ation, which predicts a much weaker suppression of heavy flavours due to the dead cone 
effect. However, coUisional energy loss is also important, as demonstrated recently by 
the Soft Colour Interaction Jet Quenching model. Based on this model we show that 
coUisional energy loss of charm quarks gives a suppression of electrons of the observed 
magnitude, but the contribution from beauty decays, which dominates at intermedi- 
ate and large pxi results in a somewhat weaker overall suppression of electrons than 
observed. 



1. Introduction 

In deep inelastic ep and pp scattering hard scattered quarks (and gluons) have 
to traverse the proton remnant, since the hard interaction takes place inside the 
proton. It cannot be expected, that interactions cease below the perturbative cut- 
off, where the coupling is large, but typical momentum transfers are small. However, 
it was founcP^, that such non-perturbative interactions may be important. As the 
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hard parton travels through the coloured background of the proton remnant, it may 
interact with it and although these soft interactions do not change the kinematics, 
they do change the colour topology. This can lead, for example, to diffractive events 
with large rapidity gaps. This circumstance is modeled by soft gluon exchanges in 
the successful phenomenological SCI (Soft Colour Interaction) model^^, but can also 
be understood theoretically in terms of multiple gluon exchange as described by the 
Wilson line in the parton density definition^. 

Following the ideas behind the SCI model it seems plausible that if soft scat- 
tering is important in DIS, it may also be important for a hard-scattered quark or 
gluon traversing a quark-gluon plasma (QGP). In the QGP many more coloured 
scattering centres are available, so that even a small momentum transfer per scat- 
tering adds up to a sizeable effect, since the interaction probability is large so that 
one parton may experience many interactions. 

2. The SCI Jet Quenching Model 

The SCI Jet Quenching modePlhas emerged from the SCI model and describes the 
rescattering of an energetic quark or gluon via soft interactions in a QGP. 

The QGP is represented by an ideal relativistic gluon gas (quarks are ne- 
glected because the plasma is initially gluon rich and the gluons come faster 
into thermal equilibrium). The gluon and energy density are accordingly given 
by Ug — 5gT^2 C(3)/27r^ and eg — 7r^f/gT''/30 respectively (gg = 16 is the de- 
generacy). The plasma expands longitudinally (Bjorken expansior3), so that the 
time dependence of the gluon density and temperature is given by Tigir) cx r^^ 
and T(r) oc t~^/^ , where r — yjt'^ — is the proper time. This rapid decrease of 
the gluon density plays an important role for the phenomenology of the model. As 
initial condition the energy density is fixed at — Ifm/c, although the plasma 
formation time can be shorter than that (r = refers to the instant of maximum 
overlap of the nuclei) . The energy density profile is proportional to the number of 
participating nucleons per transverse area npart(a^, y), since the initial energy pro- 
duction is presumably dominated by soft processes and should therefore scale with 
the number of participants. The parameter eq (5-5 GeVfm"^) gives the mean energy 
density for a central (6 = 0) collision at r = 1 fm/c, this fixes the normalisation of 
the energy density distribution for all centralities. 

The number of binary collisions and participants and the geometrical cross sec- 
tion as function of the impact parameter are calculated in the framework of a simple 
Glauber model"^ using different potentials (sharp sphere and Woods-Saxon). 

The nucleons may build up transverse momentum in soft scatterings prior 
to the hard interaction (Cronin effect). This is included by increasing the width 
(j\J^x^y,h) = -|- a ■ {Nscatix,y,b) — 1) of the Gaussian primordial k±- 

distribution of the partons by an amount a for each soft rescattering. The pa- 
rameter a — 0.5 GeV^ is obtained by fitting d-|-Au data Rhic and does not depend 
on the hard scale, since parton showers are treated explicitly. 
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The perturbative part and the hadronisation is simulated using Pythia6. 
The hard partons are produced from LO matrix elements (for heavy quarks massive 
matrix elements with nic — 1.5GeV, mb — 4.8 GeV are used) and initial and final 
state parton showers are simulated. Each of the partons is then tracked through 
the QGP and may interact with it. The interactions start at the plasma formation 
time Ti (0.2fm/c) and stop when the local energy density drops below the critical 
temperature T^. (175 MeV) or the parton leaves the overlap region. In the QGP 
a parton interacts with each plasma-gluon (with effective mass mg = 0.2 GeV) 
that is closer than a screening length i?scr = 0.3 fm (alternative parameter value 
0.5 fm) with a probability pq = 0.5 for quarks and Pg — 0.75 for gluons having 
a larger colour charge. The interactions are elastic scatterings where the squared 
momentum transfer t has a Gaussian distribution with width at = 0.5 GeV^. We do 
not use perturbative QCD matrix elements since the typical momentum transfer is 
a few hundred MeV, i.e. in the non-pcrturbative regime. Depending on the value 
of the screening length parameter i?scr being 0.3 fm or 0.5 fm, this results in a soft 
scattering cross section cr*"^' = 1.9 mb or 5.2 mb. 

After leaving the overlap region the partons are hadronised using indepen- 
dent fragmentation , since the colour connections among the partons produced 
in the perturbative phase are destroyed. A more detailed discussion of the model 
is available and we note that the parameters are the same for light and for heavy 
flavours. 

3. Results 



The results are shown for Au+Au collisions at v s^'^ = 200 GeV at mid-rapidity 

{\rj\ < 0.35). A p+p reference is simulated within the same framework and is in rea- 
sonable agreement with measurements. Gharm and beauty are simulated separately 
and added according to the production cross sections. 

A summary of the model resultJ^l for the nuclear modification factor i?AuAu for 
neutral pions in shown in Fig. [TJ With the small scattering cross section the SCI 
Jet Quenching model can account for roughly 50 % of the observed 7r° suppression 
in central collisions and comes closer to data for more peripheral events. With 
the increased cross section it is in agreement with the data for peripheral and 
very central collisions, but the increase with centrality is not quite linear as in the 
data. i?AuAu increases weakly with transverse momentum pj^, see details in our 



Fig. [2] shows a comparison of the momentum loss in a single scattering for 
different quark species. The heavier the quark is the less energy (or momentum) 
it looses. This is simply due to kinematics, since a heavy object does not take as 
much recoil as a light one. The negative part of Ap shows that low energy partons 
may gain momentum if their energy is smaller than the thermal energy of the 
surroundings (Fig. [2] shows the first scattering, so the temperature is quite high). 

In Fig.[3]we show the results for the suppression of electrons (and positrons) from 
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Fig. 1 . Nuclear modification factor for neutral pions averaged between 3 GeV and 8 GeV transverse 
momentum as function of the centrality of the collision (fraction of geometric cross section), 
Phenix data and SCI Jet Quenching Model (Woods-Saxon potential) with two different scattering 
cross sections. 
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Fig. 2. Momentum loss in the first scattering of different quarlt flavours versus their momentum. 

heavy flavour decays. While the model gives the right suppression of tt'^, with the 
same parameters the electron suppression is somewhat weaker than measured by the 
experiments. The beauty suppression is clearly weaker than the charm suppression 
as one would nai'vly expect because of the smaller energy loss. Since the electron 
spectrum at intermediate and high p± is dominated by the beauty contribution, 
the sum is relatively close to the beauty result. There is a substantial uncertainty 
on the K-factor for charm and beauty, but in this ratio only their difference is 
important. Given that pQCD calc ulat ions indicate that the charm K-factor may 
be up to twice the one for beautjff^l^ Yig. [3] show the results for K-factor ratio 
of unity or two. Fortunately, the result depends only weakly in these differences 
between the K-factors, since at large p± the sum is so much dominated by the 
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Fig. 3. Nuclear modification factor of electrons from heavy flavour decays in central collisions, 
model calculations (cr'''^' = 5.2 mb, Woods-Saxon potential) are shown for electrons from charm 
and beauty decays separately as well as the sum: dashed line for charm and beauty having the 
same K-factor, solid line for charm K-factor twice the beauty one. The model result for the nuclear 
modification factor of neutral pions is also shown for comparison. 



beauty contribution. 




P± [GoV] 

Fig. 4. of quarks that fragment into a vr" with = 5 ± 0.25 GeV (light flavours including 
gluons) or into a with px = 5 ± 0.25 GeV (heavy flavours). 

The nuclear modification factors of electrons from charm decays and neutral 
pions are similar, in spite of the smaller energy loss of charm quarks. Moreover, 
the slope of the p± spectrum of charm and the light flavour cocktail is nearly the 
same at intermediate and high p±. Fragmentation and decay affect i?AuAu; pions 
at p± — 5 GeV stem from a quark or gluon with on average 7 GeV transverse 
momentum and electrons of the same typically originate from a charm quark 
of a higher transverse momentum of 9.6 GeV, see Fig.[4l Together with the steeply 
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falling p± spectrum this can explain the similar suppression of tt*^ and electrons 
from charm decays. The beauty parent quarks lie between the other two, but the 
slope of the p± spectrum is again similar, so that the suppression of electrons from 
beauty decays is weaker. 



4. Conclusions 

In order to clarify the question, to what extent energy loss through scattering can 
help to understand the strong suppression of electrons from heavy flavour decays, we 
present results from the SCI Jet Quenching model with parameters chosen such that 
it reproduces the full vr" suppression. The main effect in the observed suppression 
of electrons can then be accounted for, but the model does not give large enough 
suppression at large p±. This is due to the weak suppression of electrons from 
beauty decays, while the suppression of charm-electrons seems sufficient. Thus, a 
pure coUisional energy loss model cannot account for the full electron suppression. 
The situation will be worse in models that include also radiative energy loss, whi ch is 
weaker for heavy flavours. The difference is here even stronger than in our modeP^. 
Since the contribution of coUisional energy loss would have to be reduced in order 
to get the correct amount of tt*^ quenching, the suppression of electr ons cannot be 
stronger than in our model (in fact a similar conclusion is drawn i It seems, 
therefore, that the conventional energy loss models miss some important part of 
physics for the heavy flavours. 
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